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ABSTRACT 

The magnetic field structure in gamma-ray burst (GRB) outflows is of great interest as it can provide 
valuable clues that might help pin down the mechanism responsible for the acceleration and coUimation 
of GRB jets. The most promising way of probing this magnetic field structure is through polarization 
measurements of the synchrotron emission from the GRB ejecta, which includes the prompt 7-ray 
emission and the emission from the reverse shock. Measuring polarization in 7-rays with current 
instruments is extremely difficult: so far there is only one claim of detection (a very high degree of 
linear polarization in GRB 021206) which despite the favorable conditions remains highly controversial. 
The emission from the reverse shock that propagates into the ejecta as it is decelerated by the ambient 
medium peaks in the optical on a time scale of tens of seconds (the so called 'optical flash') and 
dominates the optical emission up to about ten minutes after the GRB. Unfortunately, no polarization 
measurements of this optical emission have been made to date. However, after the reverse shock 
finishes crossing the shell of GRB ejecta, the shocked ejecta cools adiabatically and radiates at lower 
and lower frequencies. This emission peaks in the radio after about one day, and is called the 'radio 
flare'. We use VLA data of radio flares from GRBs to constrain the polarization of this emission. We 
find only upper limits for both linear and circular polarization. Our best limits are for GRB 991216, 
for which we find 3 a upper limits on the linear and circular polarization of 7% and 9%, respectively. 
These limits provide interesting constraints on existing GRB models. Specifically, our results are hard 
to reconcile with a predominantly ordered toroidal magnetic field in the GRB outflow together with 
a 'structured' jet, where the energy per solid angle drops as the inverse square of the angle from the 
jet axis, that is expected in models where the outflow is Poynting flux dominated. 

Subject headings: gamma-rays: bursts — polarization — radiation mechanisms: nonthermal 



1. introduction 

The detection of linear polarization at the level of 
~ 1% — 3% in the optical aftergl ow emission of sev- 
eral gamma-ray bursts (GRBs; see iCovino et alJ 120031 
for a review) has been widely considered as a confirma- 
tion that synchrotron emission is the dominant radiation 
mechanism, at least in the afterglow stage. Synchrotron 
radiation is also believed to be the dominant emission 
mechanism in the prompt 7-ray emission and in the emis- 
sion from the reverse shock, although the observational 
support for this is not as strong as for the afterglow. 
Soon after the first de tection of li nea r jjolarization in 
the afterglow emission (ICovino et al.l 1999; 'Wiiers"199!D 
it has been realized that the temporal evolution of the 
polarization (both the degree of polarization P and its 
position angle Op) can probe the magnetic field structure 
in the emitting region, as well as the structure and the 
dyna mics of GRB jets ijGhisellini &: Lazzatilll99^ iSaril 

While the polarization properties of the after- 
glow emission have received relatively wide atten- 
tion i GhiscUini & Lazzati 1999; Granot & Konigl 200^ 
yruzinov & Waxman 1999; loka & Nakamura 2001; 
Loeb f: Perna.l99R:.Medvedev fc Loebl9_95: Rossi^et^I 
■ISarilll 



tection of a very high degree of linear polarization 
(P = 80% ± 20%) in the prompt 7 -ray emission of 
GRB 021206 (ICoburn fc Boggsl f2003^ . Although this 
detection is highly controversial ( Rutledgc & Fox 200^ 
iWigger et alJl200^ . it has dramatically raised the in- 
terest in the polarization properties of the prompt 
emission. Two main explanations have been sug- 
gested for the production of tens of percent of polar- 
ization in the prompt 7-ray emission: (i) synchrotron 
emission from an ordered magnetic field in the ejecta 
llCoburn & Boggall2?M E rano ^l20M IGranot & Konigll 



l2003t " ^^^^ov!P£^e'v fc BlandfordI t2003ir and (ii) a 
viewing angle just outside the sharp edge of a jet, 
< 6'obs < do + l/Fo (Gruzinov 1999; Waxma3 
l2003j) . where 6q and Fq are the initial half-opening an- 
gle and Lorentz factor of the jet, respectively. The 
second e xplanation can work with either synchrotron 
emission f Granot' "200^; 'Nakar. Piran & Waxman' "2001 
IWaxman 2003) or invcrsc-Compton scattering of external 
photons lE ichlcr & Levinson 2003; Lazz ati et al.ll2004'at 
[Sh aviv fc Dar..l995.) . In all the above cases, the polariza- 
tion can be a good fraction of the maximal polarization 
for synchrotron emission, Pmax ~ 60% — 70%.''' 
The second explanation requires a narrow jet, OoTq < 



2004 ISaril Il999l) . the polarization of the prompt GRB 
emission received very little attention'^ before the de- 
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^ IGranot fc Kon igl (2003) had, however, pointed out that an 



ordered magnetic field in the ejecta could produce a high degree 
of polarization, P < 60%, in the emission from the GRB ejecta, 
which includes the prompt 7-ray emission, and the emission from 
the reverse shock (the 'optical flash' and 'radio flare'). 

* For inverse Compton scattering of external photons, the local 
polarization from a given point on the image can approach 100%, 
however the averaging over the unresolved image reduces the ob- 
served polarization to values only slightly higher than those for 
synchrotron emission. 
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a few, in order to have a reasonable probability of view- 
ing the jet at an appropriate angle, 9o < 0obs ^ 6'o + l/ro- 
Since a larger value of ^nrn is inferred for most GRBs 
ijPanaitescu fc Kumaill2001albr) . this would imply little 
or no polarization in most cases, where our viewing an- 
gle is inside the jet, 6'obs < ^o- The first explanation, 
however, would imply a high polarization in all GRBs, 
if indeed they have a magnetic field that is ordered over 
angular scales > I/Fq. 

The best way to reliably measure the polarization of 
the emission from the GRB ejecta is by using the emis- 
sion from the reverse shock, which includes the optical 
flash and the radio flare. While no polarization mea- 
surements have been made so far in the optical flash, 
the polarization of the radio flare can be inferred from 
existing data. In ^ we discuss the existing candidates 
for radio flare emission in GRBs and the evidence that 
this emission is indeed from the reverse shock. In |31 
we use archival data to derive upper limits on the linear 
and circular polarization of the radio flare emission, and 
in 21 show that external propagation effects in the 
host galaxy or circumburst environment are not likely 
to cause a significant depolarization the radio emission. 
The upper limits on the polarization are contrasted with 
the predictions of different theoretical models in fjs] In ^H] 
we discuss propagation effects inside the source. We dis- 
cuss the magnetic field configuration in the GRB ejecta 
in light of our results in i}71 and give our conclusions in 

m 

2. IS THE RADIO EMISSION FROM THE REVERSE SHOCK? 

In order to draw conclusions regarding the magnetic 
field structure in the GRB ejecta from the polariza- 
tion measurements of the radio flares, it is important 
to have some confidence that the radio emission is in- 
deed from the original ejecta, whose electrons were 
heate d by the reverse s hock and then cooled adiabati- 
cally (?S ari fc PiranlI999l). The best can didates are GRBs 
990123 (Kulkar ni et al.ll999aH. 991216 l|Frail et all2000(l 
and 020405 (Bcrger et al."2003'). A radio flare was also 
reported for GRB 970828 (Djorgovski ct al. 2001), how- 
ever in this case the source was detected in the radio at 
only one epoch and at one frequency, and there is no good 
evidence that suggests that this emission arises from the 
reverse shock. Also, the signal to noise in this case is 
rather bad, and does not give meaningful constraints on 
the polarization. Therefore, we shall concentrate on the 
other three GRBs for which there is better evidence that 
the radio flare emission is indeed from the reverse shock, 
and for which the quality of the data allows us to place 
interesting limits on the polarization. 

In GRB 990123 the early radio emission (up to a day 
or two) has been found to agree very w ell with the ex- 
pectat i ons for the reverse shock emission ( [Kulkar ni et al.l 
l1999at INa,ka,r fc Pir^l2nnl ISari fc P iran'l 999*). espe- 
cially when also taking int o account t he prompt optical 
emission from this GRB (jAkerlof et al.lll999|) which is 
also successfully explained as emission from the reverse 
shock. Furthermore, it is hard to explain the early radio 
emission as arising from the forward shock. Therefore, 
in GRB 990123 we have good reason to believe that the 
8.46 GHz radio measurement at t = 1.25 days that is 
used in derive limits on the polarization was indeed 
dominated by emission from the reverse shock. 



In GRB 991216 the radio emission in the first few days 
is inconsistent with the forward shock emissio n that is re- 
sponsiblc for the optical and X-ray emission ijFrail et alJ 
[2000), and a diff erent emission component is required. 
IfVail et al.i l)2000 >) suggested either reverse shock emis- 
sion, which naturally explains the early radio emission, or 
a two component jet model, where the early radio emis- 
sion is from a spherical component with Eiso ~ 10^* erg 
running into a very low external density, n ^ 10~* cm~'^. 
Such a low density is inconsisten t with bro ad band after- 
glow fits to GRB 991216 (^Panaitescu fc Kuma r 2001a ij 
which give a density higher by ~ 4 — 5 orders of magni- 
tude. Also, the energy in the spherical component is very 
high, ~ 10^"^ — lO^'' erg, which is 2 — 3 orders of mag- 
nitude higher than the values inferred for other GRBs 
jFra il et al . 2001; Panaitcscu fc Kumar 2001b). Thus, 
the reverse shock is probably the best candidate for the 
early radio emission in GRB 991216, although it is hard 
to conclusively rule out other explanations. 

In GRB 020405 the early rapid decay and spectral 
slope in the radio, F^, oc ;;-i-2±o.4j^-o.3±o.3^ provide a rea- 
sonably good case that the early radio emission (within 
the first few days) is dominated by the reverse shock 
ijBerger et alJl2'003') . 



3. LIMITS ON THE POLARIZATION FROM RADIO FLARES 

We use archival VLA observations in order to measure 
the polarization of the radio flare emission discovered 
from GRBs to date. The VLA observations are typically 
short observations (~30 minutes) that were squeezed into 
the VLA schedule for these targets of opportunity. As 
such, there is insufficient calibration information from 
these short runs alone to allow for good polarization cal- 
ibration in the standard manner using a calibrator ob- 
served over a wide range in parallactic angle. To over- 
come this difficulty we combined the VLA observations 
of the afterglows with calibration data taken from 1-4 
days before or after the run. The instrumental leakage 
terms are stable on time-scales of months, so we are con- 
fident that an optimal calibration has been obtained. As 
a check, the polarization calibration was applied to the 
absolute fiux calibrator observed on the same day as the 
radio flare (3C286 or 3C138), and the results were found 
to be consistent with their well-known properties. 

We obtain only upper limits on the linear and circular 
polarization of three radio flares. Our results are sum- 
marized in Table 1. It is worth noting that so far only 
upper limits have been found for the radio polarization 
from GRB afterglows (Tavlor ct al. 2004, 2005). In the 
best c ase of the bright G RB 030329, the 3 a limits are 
<1% (iTavlor et a,ljl200l . 



4. PROPAGATION EFFECTS EXTERNAL TO THE SOURCE 

Propagation effects might reduce the intrinsic linear 
polarization below detectable levels. A Faraday screen 
produced by ionized gas and magnetic fields can cause 
gradients in the observed polarization angle across the 
source, leading to depolarization if the intrinsic source 
size or the resolution element of the telescope is large 
compared to the gradients. The rotation measures (RM) 
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Table 1. Limits on the Polarization of Radio Flares in GRBs 



GRB 


t (days) 


tj (days) 




nc 


(MJy) 


990123 


1.25 


« 2 


< 23% 


< 32% 


242 ± 26 




1.49 




< 11% 


< 17% 


946 ± 56 


991216 


2.68 


~ 2 


< 9% 


< 15% 


634 ± 26 




1.49, 2.68 




< 7% 


< 9% 


715 ± 25 


020405 


1.19 


~ 1 - 2 


< 11% 


< 19% 


492 ± 29 



Note. — The 3 a upper limits on the linear polarization, 11^ = 
{Q'^+U'^)^/'^ / 1 , and circular polarization, lie = V/Iy of radio flares 
from GRBs, where Q, U, V and I are the Stokes parameters. In 
the third line for GRB 991216 we combine the two epochs from the 
first two lines in order to obtain a better limit on the polarization. 
These limits were derived using VLA observations at 8.46 GHz. 

can be related to the line-of-sight magnetic field, , by 



where the upper limit of integration, L, is the distance 
from the emitting source to the end of the path through 
the Faraday screen along the line of sight. To produce a 
90° rotation at 8.4 GHz requires a RM of 1200 rad m^^ 

In our Galaxy the RM can reach up to ^ 
1000 rad m~^ at very low Galactic latitudes, while at 
high Galactic latit udes it drops to a few x 10 rad m~^ 
dSimard-Normandi n. Kronberg fc Buttonlll98lD . At low 
Galactic latitudes iGlegg et al.l l|199^ find a RM differ- 
ence of up to 180 rad m~^ on scales of 1° which trans- 
lates to a linear scale of ~ 100 — 200 pc and gradient 
of ARM/AL_L - 1 radm"^ pc"^ in the RM. Since in 
our case the source size is ^ (1 — 3) x 10^^ pc, this 
would correspond to a negligible rotation measure gradi- 
ent across the image, ARM ~ (1 — 3) x 10^^ rad m^^, 
some 5 orders of magnitude too low to cause signifi- 
cant depolarization. The host galaxies of GRBs are not 
thought to be greatly different from our own Galaxy 
ijBloom. Kulkarni. fc Diorgovskill200^ . 

Molecular clouds in our Galaxy typically exhibit vari- 
ation of ARM ~ 18 - 30 rad m"^ in t he RM over 
scales of - 2 pc (iWolleben fc Reichll200^ . This cor- 
responds to ARM/ALj_ ~ 9 — 15 rad m~^ pc""'^ and 
ARM ~ 0.1 — 0.4 rad m~^ across the image, which is 
still some 4 orders of magnitude too low to cause signif- 
icant depolarization. Individual Hii regions can produce 
enhanced RMs with dispersions of order ^ 50 rad m~^ 
on scales of ~ 0.2 pc (jGaensler et al.ll200l|) . correspond- 
ing to ARM/ALx ~ 250 rad m-^ pc"! and ARM ~ 
2.5 — 7.5 rad m^^ across the image, which is still some 
2-3 orders of magnitude too low. 

Now we consider propagation effects in the immediate 
environment of the GRB, which was shaped by its pro- 
genitor. If the latter is a massive star, the immediate 
environment is the pre-explosion stellar wind which for 
a constant mass loss rate and wind velocity would give 
a density profile next oc R~^, at radii R that are smaller 
than that of the wind termination shock. Such a density 
profile would imply a small deceleration radius i?dcc and 
a very high density at that radius, which would in turn 
imply fast cooling of the reverse shock electrons and no 
detectable radio flare emission at t ~ 1 day. Keeping this 
in mind, we shall assume a uniform external medium. 



For GRB 020405. ISeTger et all l|2003|) model the opti- 
cal. X-ray and radio light curves with a uniform medium 
of density next ~ 0.05 cm~^. Afterglow models applied 
to GRB 020405 predict an intrinsic diameter of the radio 
emitting region at 1.2 days after the burst of ^ 10^^ cm. 
Values of the external density inferred from broad band 
afterglow modeling ca n be as high as n < 30 cm~^ 
ijPanaitescu fc Kumad l2001albl: lYost et alJ 120031) . For 
GRB 990123 a very low density of next ~ 10"^ cm^^ 
is inferred, while for GRB 991216 the inferred density is 
more typical, next 5 cm^"^. 

Magnetic fields within the circumburst medium are not 
expected to cause a large depolarization, unless there is 
a clump with a high density contrast along our line of 
sight. A clump much smaller than the image size would 
have a small effect. A clump larger or comparable to the 
image size, Lim, would produce a change ARM in the 
rotation measure across the image which is roughly inde- 
pendent of the clump size, L, for a given density contrast 
C. This can be seen as follows. We have n = Cnext oc C 
and assuming flux freezing and isotropic compression, 
B cx n2/3 oc so that RM oc nBL oc C^/^L, 

and the change in the RM across the image is roughly 
ARM ~ RM X Lini/L oc C^/^, where the dependence on 
L cancels out. Furthermore, the gradient in the rota- 
tion measure across the image would be comparable to 
the total rotation measure of a clump with a size similar 
to that of the image, ARM ~ RM(L = Lim), causing a 
relative rotation of the polarization position angle of 

^ < 45 (^^) (^) (^^] f^V' 
90° ■ Vlcm-3/ VlOyuGy/ VlQi^cmy VlOOy 

(2) 

where next and i?ext are the number density and magnetic 
field of the un-clumped external medium, and R± is the 
apparent radius of the image on the plane of the sky. 
Thus a clump with a density contrast C > 10^ and a 
size > 10^^ cm is required in order to induce significant 
depolarization. We also note that if such a clump lies 
within a distance of < 1 pc from the progenitor star it 
would produce a detectable bump in the afterglow light 
curve, which was not observed for GRBs 990123, 991216 
or 020405. Thus, such a clump along our hne of sight 
seems unlikely. It is much less likely that such a clump 
would be along our line of sight in all the three GRBs 
that we analyzed. 

Instead of a single clump, one might envision many 
small clumps with a covering factor of / ~ 1 (in order to 
cover most of the image) . The typical clump size would 
have to be Ln = L/(10^^ cm) < 1 in order to avoid 
producing detectable bumps in the afterglow light curve. 
Since even for C2 = C/100 1 we need Ln > 1 in order 
for a single clump to produce significant depolarization 
(see Eq. El, many small clumps are required along a ran- 
dom line of sight. If their rotation of the position angle 

would always be in the same direction, the required num- 

5/3 1 

ber of clumps would be Nq > C2 L^j . However, since 
the direction of rotation of the position angle is expected 
to be random between different clumps, the total rota- 
tion would on average be zero with a mean r.m.s value 

1 /2 

of Nq' times that of a single clump. Thus a strong de- 
polarization would require > Nq clumps along a random 
line of sight. Thus Nd{R)L^/R^ > N§ - C^'^^^^L^^ 
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where Nci{R) is the total number of clumps at radii 
^ R. Thus the total mass in the clumps would be 

~ O.lM0(next/l cm~'^)C2~^^'^ij"y^i?Jg while there volume 

filling factor would be ~ O.IC^^^^^ L^j^ R^^ . This means 

that the clumps would hold ^ 10C2'^^^L^jRYg times 
more mass than that of the material between the clumps. 
Such an extreme clumping seems highly unlikely. 

5. THE RESULTING CONSTRAINTS ON THEORETICAL MODELS 

The radio flare emission typically peaks on a time scale 
of '-^ 1 day after the GRB. This is usually similar to the 
jet break time, tj, in the afterglow light curve. At this 
time the Lorentz factor is typically T ^ 10, which is much 
smaller than its initial value, Fq > 100. Thus a good 
part of the jet (or all the jet for t > tj) is visible, and the 
observed polarization is an effective average value over 
this observed region. The optical flash emission peaks 
at the deceleration time idee when the Lorentz factor is 
close to its initial value (unless the reverse shock is highly 
relativistic), and should thus have a polarization close to 
that of the prompt 7-ray emission. It could therefore 
provide information that is complimentary to that from 
the radio flare. 

The jet break time, tj, plays an important role in 
the polarization light curve for most theoretical mod- 
els. Therefore, it is important to know its value when 
comparing between theory and observations. In GRB 
990123 ther e is a jet break in the optical light curve at 
tj ~ 2 days l|Kulkarni et alJll999hl: iPanaitescu fc Kumail 
l2001ajl . A similar jet break time, tj « 2 days, was found 
in th e o ptical light curve of GRB 991216 (HalDcrn et al^ 
I2OOO: Pa naitescu fc Ku mar 2001a) with a rather large 
uncertainty on this value. In GRB 020405 the jet break 
is not readily apparent in the data.^ iBerger et al.l 1)20031) 
deduce tj « 0.95 days from the fit to the broad band 
afterglow light curve, while Price et al. (2003) infer tj — 
1.67 ± 0.52 days. Thus we use a value of ~ 1 — 2 days 
for GRB 020405. In all cases the polarization limits are 
at times either similar to or slightly before the jet break 
time, t "^tj. 

If the observed radio frequency is below the self ab- 
sorption frequency, u < J^sa, this can significantly re- 
duce the polarizatio n, for any magnetic field configura- 
tion. iNakar fc Pira n (2004) found that the peak of the 
radio flare emission is typically due to the passage of the 
self absorption frequency across the observed band. This 
would imply i> < Vgn and a suppression of the polariza- 
tion during the rise to the peak, and v > Vsa. (no sup- 
pression of the polarization) during the decay after the 
peak. This can reduce the polarization before or around 
the time of the peak in the radio flare emission. In GRBs 
991216 and 020405 the radio flux is decaying at the time 
of the measurement of the radio flare, and direct mea- 
surements of the spectral slope also support an optically 
thin spectrum (y > i^sa), so that no suppression of the 

^ For example, IMasetti et al.l 120031) find no evidence for a jet 
break in the fight curve up to ~ 10 days. Since the most severe 
constraints from our polarization measurements are for an ordered 
toroidal magnetic field in the ejecta, in which case a larger value 
for the jet break time {tj > 10 days) would imply even stricter 
constraints on the model, we adopt a conservative approach and 
use the lower va l ues of tj that were inferred bv iBerger et alj 120031) 
and lPrice eTall <2003) . 



polarization is expected due to self absorption. For GRB 
990123 the point used to measure the polarization seems 
to be near the peak of the radio flare and there is no 
good measurement of the spectral slope near this time. 
Thus the polarization might be somewhat suppressed if 
V Vsa., but probably not significantly suppressed since 
the relatively high flux and the proximity to the peak of 
the radio flare suggest that v is not much smaller than 
and the two are at most comparable. 

5.1. The Ejecta Dynamics after the Reverse Shock 

Since the radio flare emission typically peaks at t < tj , 
the possible lateral spreading of the jet at i > tj can 
be neglected, and the jet dynamics can be reasonably 
approximated as being part of a spherical flow. After 
the passage of the reverse shock, the shocked external 
medium approaches the Blandfo rd fc McKc c (1976, here- 
after BM) self similar solution. However, the shocked 
ejecta has a significantly higher density than that given 
by the BM solution, and thus its rest mass density be- 
comes comparable to its internal energy density (i.e. it 
becomes "cold" ) much earlier on. For a mildly relativistic 
reverse shock this happens soon after the reverse shock 
crosses the shell, while for a relativistic reverse shock 
the shocked shell is first reasonably described by the BM 
solution, while it is still "hot" , but deviates from that so- 
lution once it becomes "cold" ftKo bav ashi_fc Sari 2000). 

For the BM solution the value of the self similar vari- 
able X for a fixed fiuid element, that is appropriate for the 
original ejecta, evolves with radius R as x= {R/Ro)'^~ '' 
for an external density pcxt o^r~^ l)Granot fc S ari 20021. 
where Rq is the radius where it crossed the shock, which 
for the original shell of ejecta is given by the deceler- 
ation radius, i?o ~ -Rdec- Also, 7 = Fx"^^^ where 
F oc i?-(3-fc)/2 ^j^g Lorentz factor of the fluid just be- 
hind the shock, so that 7 oc i?~('^~2fc)/2 ^j^g original 
ejecta. More generally one can assume some power law 
dependence of the Lorentz factor on radius, 7 oc i?^^, 
where the power law index g can deviate from its value 
for the BM solution, g = 7/2 — fc. This implies t oc 
i?/72 oc 7-(29+i)/9 OC, 7 OC tsU'^a+i) and R oc ti/(2s+i). 

Once the shell becomes "cold" (i.e. its rest energy 
exceeds its internal ener gy) its dynamics would deviate 
from the BM solution l|Kobavashi fc Saril 1200(11) . how- 
ever, the power law index g is still bounded by the value 
just behind the forward shock, g — (3 — fc)/2 (since 
the the ejecta shell is lagging behind the shocked ex- 
ternal medium) and by the value for the BM solution, 
g = 7/2 — fc (since there is a larger inertia that resists 
the deceleration which is driven by a similar pressure). 
Thus, (3-fc)/2 < 5 < 7/2-fc or 3/2 <g < 7/2 for fc = 
(a uniform density external medium).® 

At any given observed time t > tdcc the observed 
emission from the reverse shock is from a somewhat 
smaller radius and a slightly smaller Lorentz factor com- 
pared to the emission from the forward shock. This 
can be seen as follows. Let quantities normalized to 

^ A stellar wind environment is not so relevant for the radio flare 
since in this case the reverse shock electrons are are fast cooling 
(i.e. cool significantly due to radiative losses within the dynamical 
time) and therefore the observed flux rapidly decays after the pas- 
sage of the reverse shock, and no detectable radi o flare is expected 
IChevalier fc Lill200Cl: iKobavashi fc Zhanell2003) . 
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Fig. 1. — Our 3 a upper limits for the linear polarization of the 
radio flare emission overlaid on the theoretical polarization light 
curves for a toroidal magnetic field in the GRB ejecta. The error 
bars represent the uncertainty in the determination of the jet break 
time tj from the optical afterglow light curve. The to p tw o panels 
are for a uniform jet, and are calculated in AppendixKl The dif- 
ferent lines, from top to bottom are for dahs/^O = 0.9, 0.8, 0.1. 
The main figure is for a = — dlog i^^/dlog = — 1/3 and Pmax = 
(a + l)/{o -I- 5/3) = 1/2, while the inset is for a = 3/4 and 
^max = 21/29 fa 0.724. In the top panel the Lorentz factor of the 
ejecta is assumed to remain equal to that of the freshly shocked 
fluid just behind the forward shock (FS), while the middle panel 
it is assumed to follow the Blandford & McKcc (1976) (BM) self 
similar solution. The bottom panel is for a 'structured' jet where 
the energy per soli d angle drops as outside of some small core 
angle (taken from iLazzati et al.ll2004bh . In this case P(t/tj) is 
practically independent of d^^,^. 



their values at t^cc be denoted by a twiddle. From 
the equality of the observed time we have i?Rs — 
^(4-fe)/(2,+i) ^ ^(2,.s+i)/(29Hs+i) ^heregps = (3-fc)/2 
for the forward shock and g — g^is for the reverse 
shock. Since R oc ^1/(29+1) this implies that the ratio of 
the Lorentz factors for a fixed observed time is 7/r = 

n-(9RS-9Fs)/(2SRS + l) _ f-(SRS-9Fs)/[(2SRS + l)(29FS + l)] _ 

f-[3-(3-fe)/2]/[(4-fe)(23+i)]_ For fc = this givcs 7/r = 
f-(9-3/2)/[4(2g+i)] .^y]-^g;^-g ^hc exponent ranges between 
— 1/16 and for the possible range of 3/2 < g < 7/2, 
which implies 7/r « iJ The somewhat smaller radius 
of the ejecta shell and the fact that it is typically cold 
by the time of the radio flare which is typically aX t < tj 
suggest that a possible lateral spreading of the jet is prob- 
ably not important for the radio flare emission, and thus 
it shall be neglected when calculating the polarization. 

5.2. The Implications for Different Theoretical Models 

Since typically t^cc ~ 10 — 100 s and the radio flare emission 
is observed at t 10^ s, this implies 1 < r/7 < 1.5 — 1.8. 



If a high polarization in the prompt GRB is caused by 
a viewing angle < ^obs ^ ^0 + 1/ro, then the polariza- 
tion of the optical flash around its peak would be similar 
to that of the prompt 7-ray emission. This is since at 
this time (i ~ idee) T ^ Tq and Qj « (i-e- very little 
lateral expansion could have taken place before the de- 
celeration time). However, at the time of the radio flare 
(t ^ \ day) r ^ To and even a modest lateral expansion 
would give ro(0j —^0) ^ 1, so that the jet would occupy 
our line of sight. This could signiflcantly reduce the po- 
larization for a magnetic field that is random within the 
plane of the shock, as is expected to be produced in rela- 
tivistic collisionless shocks by the two stream instability 
(jMedvedev &: Loebl fl999). Furthermore, for a shock pro- 
duced magnetic field the polarization of the radio flare 
should be similar to that of the optical emission from the 
forward shock, which is observed to be P < 3%.^ 

Such a low polarization suggests that the magnetic 
field is not random in 2D, fully within the plane the 
shock, but is instead more isotropic, and has a compa- 
rable c omponent in the direction normal to the shock 
(jGrano t & Konigl 2003). Constraints similar to our best 
case, GRB 991216, for a larger number of GRB, could 
suggest a similar conclusion for the magnetic field in the 
ejecta.^ 

Our upper limits on the polarization of the radio flare 
emission can put much tighter constraints on models 
where there is an ordered magnetic field in the GRB 
ejecta. If there is a magnetic field in the ejecta that 
is ordered over patches of angular scale Qb ^ I/I'd 
then typically P ~ Pmax during the prompt 7-ray emis- 
sion and near the peak of the optical flash, while dur- 
ing the radio flare the polarization can be reduced due 
to averaging over N ~ (i9r)~^ incoherent p atches, 
P Pmax X min[l,76'B] l|Granot fc Konigll 12003). This 
implies 6*^ < Piim/7Pmax, where Pum is our upper limit 
on the linear polarization. During the radio flare 7 ~ 10 
so that for GRB 991216 we have Pum = 0.07 and 
Ob < 1.4 X 10-2(7/10)-i(Pmax/0.5)-i rad. In partic- 
ular, if the magnetic fleld is roughly uniform over the 
whole jet (Ob ^ 9o) then this would imply P ^ Pmax 
ijGranot fc Kon igl 2003), i.e. tens of percent of polar- 
ization, which is definitely inconsistent with our upper 
limits. 

For a toroidal magnetic field in the ejecta, we 
show our upper limits superimposed on the theoret- 
ical polarization light curves in Fig ^ both for a 
uniform jet and for a 'structured' jet. The lat- 
ter is expected in models where the GRB outflow is 
Poyn ting flux d ominated fLvutikov fc Blandford' '200^ 
Lvuti kov. Parle y & Blandford 2003). The polarization 
light curve for a structured jet is taken from Lazzat i et all 
(j2004h) . The polarization light curves for a uniform jet 
are derived in Appendix^ For a uniform jet we consider 
the two limiting cases for the evolution of the Lorentz 
factor of the ejecta after the deceleration time (see i )5.1l 

^ In order to detect the afterglow emission and measure its po- 
larization one must flrst detect the prompt 7-ray emission, which 
requires So^s ^ ^0 I/Tq. Since typically Fo^o 3> 1, most lines of 
sight would be inside the jet, (?obs < ^0- 

^ A line of sight sufficiently close to the jet axis could produce a 
very low polarization. However, such viewing angles correspond to 
a small solid angle, so that the probability for such lines of sight 
in all of the GRI3s in a reasonably sized sample is very small. 
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for details): (i) the same as the forward shock Lorentz 
factor, and (ii) folfowing the iBlandford fc McKej ^197(t) 
self similar solution. In case (ii) the Lorentz factor of 
the ejecta 7 decreases slightly faster with the observed 
time t (cx i"^/!*^ instead of t^^^^). This shghtly 'con- 
tracts' the polarization light curve along the time axis. 
Also, while ^Oj{tj) = 1 in case (i), we have j9j{tj) — 

(ijAdcc);'/'^ ~ 0.65 where typically i^Adoc ~ 10^ 
which shifts the polarization light curves to earlier times 
by a factor of (tj/tdoc)^^^ 2.7. This can be readily seen 
by comparing the upper and middle panels in Fig. ^ 

Since P^ax = (a + l)/(a + 5/3) (Granot 2003), a higher 
values of a produces a higher degree of polarization. Nev- 
ertheless, even their lowest values for optically thin syn- 
chrotron emission, a = —1/3 and Pmax = 1/2, still pro- 
duce a fairly high degree of polarization. 

For a uniform jet, P{t ^ tj) significantly increases with 
Qohs/Oa, and goes to zero at 6'obs — 0. Thus, our upper 
limits on the polarization put upper limits on ^obs/^o- 
These limits also depend on the dynamical model for the 
GRB ejecta. For GRB 991216 we obtain Oohs/Oa < 0.4 
and 0.55 for cases (i) and (ii), respectively. 

The model that is most severely constrained by our up- 
per limits on the polarization is a toroidal magnetic field 
together with a structured jet (see lower panel of Fig. 0). 
In this case all of our upper limits are significantly below 
the predictions of this model (by a factor of > 6 for GRB 
991216). Thus a predominantly toroidal magnetic field 
in the GRB ejecta together with a structured jet is hard 
to reconcile with our upper limits. 

6. PROPAGATION EFFECTS AT THE SOURCE 

The effects of propagation in plasma inside the 
source on the synchrotron emission in G RBs has 
been recently considered iMatsumiva fcTo ka 200^ 
ISas-iv. Wa,yma,n T>oeHl20(il hereafter SLW04). In the 
early emission while the reverse shock is still going on, 
these effects can suppress linear polarization and produce 
up to tens of percent of circular polarization below the 
self absorption frequency Vsn, for a magnetic field that 
is ordered on large scales (SLW04). We have obtained 
upper limits both on the linear polarization (11^) and on 
the circular polarization {He) for the radio flare emis- 
sion. The best limits are for GRB 991216: Ul < 7% and 
He < 9% (3 a). Our upper limits directly constrain such 
plasma propagation effects. 

There are some arguments which suggest that plasma 
propagation effect should be relatively small and sub- 
dominant in the radio fiare emission. First, the suppres- 
sion of linear polarization and prevalence of circular po- 
larization is much larger when the reverse shock electrons 
are fast cooling, i.e. most electrons cool significantly on 
a time scale smaller than the time it takes the reverse 
shock to cross the ejecta shell (SWL04), and is less sig- 
nificant when they are slow cooling (Matsuniiva & loka 
l2003j) . As mentioned in footnote |H| slow cooling is re- 
quired in order to have a detectable radio flare emission, 
and is therefore the relevant case for us. Second, the 
effects of propagation in plasma decrease with time at 
a fixed observed frequency v (fMatsumiva & loka 2003) 
and are significantly smaller at the time of the radio flare 
(t ^ 1 day) compared to the time when the reverse shock 
is still going on (^dcc ^ lO'^ s). Furthermore, 11^ < 1% 



is expected at v > v^^ (jMatsumiva fc Iokall2003(l . which 
seems to be the case for the radio flare emission that we 
consider in this paper. 

Another related effect that might cause depolarization 
(SWL04) is a different amount of Faraday rotation for 
different frequencies at a given point on the image (which 
would be hard to resolve because of the finite instrumen- 
tal spectral resolution) or at the same observed frequency 
between different points on the image (because of the 
different path lengths through the emitting plasma and 
different Doppler factor corresponding to different fre- 
quencies in the local rest frame of the emitting plasma). 
The lS.v/v which gives a Faraday rotation of 27r for a 
cold plasma is — ■jT^mlc^^v')'^ /{e^neB'W) where 

W' is the width of the emitting region in the local frame 
(SWL04). For a relativistic plasma with a typical elec- 
tron random Lorentz factor of 7e this expression should 
be multiplied by ~ (ln7e)/7g (Matsumiya fc loka 2003; 
SLW04) so that Faraday rotation and the resulting depo- 
larization are suppressed. SWL04 consider fiducial pa- 
rameters that result in fast cooling, and obtain that a 
significant depolarization is possible up to observed fre- 
quencies of ~ 10^^ Hz at t < idee- Since we consider 
v ^ 10^" Hz and Faraday rotation scales as we need 
to suppress Faraday rotation by some 10 orders of mag- 
nitude compared to their estimate. 

Since a detectable radio fiare requires slow cooling of 
the reverse shock, this would imply a smaller density 
and magnetic field in the local frame, which would re- 
duce the Faraday rotation at i < tdec- Furthermore, 
SLW04 used t^cc = 10 s while the radio fiare occurs at 
t ~ 10^ s. We have UeW oc R''^ and for the BM so- 
lution W oc i?/7 oc and B' oc l/RW oc i?-"/^ 
so that HeB'W oc _R~^^/^ oc t"^^/^^, which reduces the 
Faraday rotation by ^ 4 orders of magnitude between 
idoc and the time of the radio fiare. Also, at the time of 
the radio fiare 7 ^ 5 — 10 compared to ~ 10^-^ at idee so 
that the same observed frequency corresponds to a higher 
frequency in the local rest frame of the emitting plasma 
{v « 7J^')- This reduces Faraday rotation by ~ 3 — 4 
orders of magnitude. Finally, the reverse shock electrons 
must still be relativistic during the radio flare in order 
to emit the observed synchrotron radiation, and this will 
suppress the Faraday rotation by a factor of ~ (ln7e)/7e 
(SWL04). Altogether, it seems that over a large part of 
the relevant parameter space there will not be significant 
depolarization of the radio fiare emission due to differ- 
ent amounts of Faraday rotation at different points on 
the image, although it is hard to rule out some degree of 
depolarization for certain regions of parameter space. 

7. DISCUSSION 

If the magnetic field is carried out from the source, 
then at large radii it would be almost completely in the 
plane perpendicular to the radial direction, since the ra- 
dial component of the magnetic field scales as while 
the two transverse components scale as R~^ if the mag- 
netic field is initially tangled on small scales. If the flow 
and the magnetic fleld configuration are axially symmet- 
ric then B^ oc R~^ while Bg oc R~^ so that B^ will 
dominate at large radii and the magnetic field will be 
toroidal. If the ratio of electromagnetic to kinetic energy 
is cr < 1, then the magnetic field is sub-dominant dynam- 
ically and the plasma can in principle keep it tangled on 
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small scales within the plane normal to the radial direc- 
tion, if it is initially tangled on such small scales, and 
its exact configuration is not obvious from theoretical 
considerations. For example, it could be ordered on an- 
gular scales 9b Oj where 9j is the half-opening angle 
of the jet. On the other hand it might be ordered over 
the whole jet, possibly in a toroidal configuration. If, on 
the other hand, the GRB outflow is Poynting flux dom- 
inated (i.e. (T 3> 1) then the magnetic field is dynam- 
ically dominant and can move the plasma around and 
arr ange itself i n a toroidal configuration over the whole 
jet l|Lvutikov fc Blandford 2003, 2004). Thus an ordered 
toroidal magnetic field is expected in the latter type of 
model. 

These are the magnetic field configurations that are 
expected before the prompt gamma-ray emission. For 
cr < 1 the gamma-ray emission is likely due to inter- 
nal shocks, which may cause a shock produced magnetic 
field that is expected to be random within the plane of 
the shock, so that it can add a random component to an 
initially ordered component. A similar effect can also oc- 
cur in the reverse shock. The random component (-Bind) 
can reach values close to equipartition, so that it might 
be dominant over the ordered component (Boid) if the 
latter is well below equipartition (ct ^ 1). However, for 
(T ~ 1 with an initially ordered magnetic field, the shock 
produced random component can at most be comparable 
to the ordered component. 

In Poynting flux dominated models, where initially 
(7^1, there would be no reverse shock and therefore 
no radio flare if a remains ^ 1 after the prompt gamma- 
ray emission. Thus, in order to get a radio flare we need 
cr < 1. This might still be possible in such a model 
if after the dissipation of electromagnetic energy (mag- 
netic reconnection) that causes the gamma-ray emission, 
the value of a decreases to ~ 1. Furthermore, the dissi- 
pation of magnetic energy (magnetic reconnection) that 
gives rise to the prompt gamma-ray emission in this type 
of models can also change the local configuration of the 
magnetic field in the emitting region, so that it might not 
be perfectly toroidal, and could also have a random com- 
ponent that could potentially be comparable in strength 
or possibly even exceed the ordered component. 

A combination of an ordered and a random mag- 
netic field component could decrease the resulting de- 
gree of polarization, compared to that for a purely or- 
dered magnetic field , by a factor of ~ y/jl + rj) where 
V «^ (-B^rd)/(^r„d) HCTranot fc Konigll 12001 . Thus we 
obtain r] < 0.2 for GRB 991216, where our upper limit 
on the polarization is > 6 times smaller than the pre- 
dicted value for a purely toroidal magnetic field. In other 
words, a sufficiently sub-dominant ordered toroidal mag- 
netic field together with a larger random magnetic field 
component is required. This requirement is not trivial 
when starting from an ordered toroidal magnetic field 
close to the equipartition value, since the shock produced 
random field component would not exceed equipartition 
and thus would typically be at most comparable to the 
ordered component (77 > 0.5). 

8. CONCLUSIONS 

We have derived upper limits on the linear and cir- 
cular polarization of the radio flare emission discovered 
from GRBs to date. Our results are summarized in Table 



1. There is a reasonably good case that the radio flare 
emission in GRBs 990123, 991216, and 020405 indeed 
arises from the original ejecta that was shocked by the 
reverse shock and then cooled adiabatically (as discussed 
in Sj21l- This emission is also most likely predominantly 
synchrotron radiation. Therefore, our upper limits on 
the polarization can be used to constrain the magnetic 
field structure in the ejecta of these GRBs. 

Propagation effects outside of the emitting region 
might decrease the measured polarization below the 
value of the intrinsic polarization (see discussion in Sj^J. 
We have demonstrated that the relatively small source 
size, < 10^^ cm, at the time of the radio flare makes 
it very unlikely that gradients in the rotation measure 
across the image due to the magnetic field in the ISM of 
the host galaxy would cause significant depolarization. 
Depolarization due to a possible propagation in a molecu- 
lar cloud or in the immediate circumburst medium might 
cause larger gradients in the rotation measure across the 
image, but still require extreme conditions in order to 
cause significant depolarization. 

Propagation effects in the plasma inside the source 
might cause a different amount of Faraday rotation at 
different points on the image which may cause depolar- 
ization since the image is not resolved. For a fixed ob- 
served frequency, such a depolarization is much smaller 
during the radio fiare compared to t < idee (see discus- 
sion in Since Faraday rotation scales as it might 
still be important in the radio for some parameter val- 
ues, although it should not be very important for a large 
part of the relevant parameter space. Thus a more thor- 
ough investigation of this effect and its magnitude over 
the relevant parameter space is called for. Keeping this 
caveat in mind, we continue and examine the implica- 
tions of our upper limits on the polarization under the 
assumption that there was no significant depolarization. 

We have compared our upper limits to the predictions 
of different theoretical models. Models where the mag- 
netic field is produced in the shock would in most cases 
produce a polarization below our upper limits. Further- 
more, they are expected to produce a polarization similar 
to that of the afterglow emission, if indeed the dominant 
magnetic field in the afterglow is shock produced and 
if the magnetic field configuration behind the afterglow 
shock is similar to that behind the reverse shock. For 
lines of sight near the edge of the jet the expected polar- 
ization might in some cases exceed our upper limits, if the 
magnetic field behind the shock is maximally anisotropic 
(i.e. random within the plane of the shock or ordered 
in the direction normal to the shock). This might sug- 
gest a more isotropic magnetic field configuration behind 
the shock if comparable or better upper limits will be 
found in a larger sample of GRBs, similar to the current 
situation with the afterglow emission where linear po- 
larization has been measured in several GRBs and was 
found to be 11^ < 3% in all cases, perhaps with one 
exception - GRB 020405 where a sharp spike in the po- 
larization (P = 9.9 ± 1.3% at t — 1.3 days) was reported 
(jBersier et al.ll200^ .^° We note that an almost simul- 

This optical polarization spike occurred at a similar time to 
the measurement of the radio flare emission from the same GRB 
which we used in order to derive the upper limit on its polarization 
(see Table 1). This is probably a coincidence. 
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taneous polarization measurement by a different group 
ijMasetti et al.ll2C)?)^ resulted in a significantly lower po- 
larization, P = 1.5 ± 0.4%. 

Our upper limits on the polarization put much stronger 
constraints on models in which there is an ordered mag- 
netic field in the ejecta. A magnetic field that is roughly 
uniform across the whole jet would produce tens of per- 
cent of polarization, which is inconsistent with our upper 
limits. If the magnetic field is ordered over patches of an- 
gular scale 6b, which are mutually incoherent, then our 
upper limits on the polarization put upper limits on Ob- 
The tightest constraint is for GRB 9901216, for which 
we find Ob < 10^^ rad (see C!^ . 

The polarization light curves for a toroidal magnetic 
field in the ejecta, together with our upper limits, are 
shown in Fig. ^ For a uniform jet, our upper limits 
on the polarization constrain our viewing angle toward 
GRB 991216 to be 6'obs/6'o < 0.4 - 0.55 where this range 



roughly covers the uncertainty in the dynamics of the 
ejecta at t > tdcc (as discussed in ii5.2|) . These values are 
for a conservative value of the spectral slope, a = —1/3 
which corresponds to Pmax — 1/2. For a structured jet 
with a toroidal magnetic field, P Pmax is expected at 
t ^ tj (which applies for all the upper limits given in Ta- 
ble 1). Therefore, this model appears to be inconsistent 
with our upper limits on the polarization. 
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APPENDIX 

POLARIZATION OF A UNIFORM JET WITH A TOROIDAL MAGNETIC FIELD 

Here we derive the polarization of a uniform jet with an ordered toroidal magnetic field, by generalizing the results 
of IGrana3 l)2003|) for a uniform transverse magnetic field. The emission is assumed to arise from a section of a thin 
spherical shell moving radially outward with a bulk Lorentz factor 7^1, that lies within a cone of half opening 
angle Oj, which represents the jet. For simplicity, the emission is integrated over the jet at a fixed radius and the the 
differences in photon arrival time from different angles 9 from the line of sight are ignored. An integration over the 
equal arrival time surface of photons to the observer might introduce small qu antitative diffe rences, but the results 
should be qualitatively similar, as we verified by comparing our results to those of lLazzati et al.l (j2004b). The emission 
in the local rest frame of the shell (where quantities are denoted by a prime) is taken to be uniform across the jet 
(hence a uniform jet), and depends only on the angle x' between the direction of the emitted radiation, n , and the 
local direction of the magnetic field, B'. The polarization position angle at a given point on the jet makes an angle of 
&P,B = + arctan{[(l — y)/(l + y)] cot(f>} from the local direction of the magnetic field, B fGran ot fc Konigill2003(l . 
where y = (7^)^, and (p is the angle between B and the direction from the line of sight to that point on the jet. 

The Stokes Parameters are given by (J7,Q)//Pmax = / dnii,{sm29p, cos29p)/ J dill^, where 9p is measured from 
some fixed direction, which for convenience we choose to be the direction from the jet symmetry axis to the line of 
sight. We have dfi oc dipdy where tp is the azimuthal angle around the line of sight measured from the direction 
between the jet axis and the line of sight. Let us use the notations q = 6'obs/^07 Vj = il^j) 
a = 9/9obs = q-Hylyjf'^. We have h = I'Av/v'f with J^, cx (z/')-"(sinx')', v/v' « 27/(1 



sm' 



x' = l-(n'-B') 



1 



y 



cos 



1 



y 



y± = (1 ± <?) yj, and 
I- y) and 

cos 



(1 + yY [l + a? + 2a cos (^) 



(Al) 



where n' is the direction in the local frame of a photon that reaches the observer. We also find that 

9p ^ (f — arctan 





-,\ 


sin (/J 




+ yj 


(a -I- cos ip) _ 



(A2) 

From symmetry considerations U = Q and P = \Q\/I. The direction of the polarization vector on the plane of the sky 
is along the line connecting the jet symmetry axis and the line of sight. The degree of polarization is given by 
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where Q{x) is the Heaviside step function, cos 'Si = [(1 

producing the results shown in Fig. ^we use e — 1 + a ljGra,noHl200^ \ 

In order to produce polarization light curves a simple model of j(t) is added, where t is the observed time. For a 
jet with no lateral spreading going into a uniform density medium 7 oc where ^ = 3/8 + {g — 3/2)/[A{2g + 1)] and 
7 oc (see ^5.1|l . We identify the jet break in the optical light curve with r9j = 1 where F is the Lorentz factor 
just behind the forward shock. Thus 



1/2 a 

y/ =l^j 



[g-(3-fe)/2]/[(4-fe)(2g+l)] 



-g/(2g+l) 



(A4) 
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As discussed in >i5.1l we have (3 — k)/2 < g < 7/2 — k. For the lower hmit 5 = (3 — k)/2 which corresponds to 7 = F 
we have = {t/tj)-'~^-''^l^'^-^\ For the upper Umit 5 = 7/2 - A: we have = (tj/idoc)^/'^''^"''^! (V^j)"^^"^''^^'^^'^"''^' • 
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